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ON A SOLUTION OF THE NONLINEAR DIFFERENTIAL EQUATION FOR TRANSONIC FLOW 
PAST A WAVE-SHAPED WALL ' 


By CARL KAPLAN 


SUMMARY 


The Prandtl-Busemann small-perturbation method is utilized 
to obtain the flow of a compressible fluid past an infinitely 
long wave-shaped wall. When the essential assumption for 
transonic flow (that all Mach numbers in the region of flow are 
nearly unity) is introduced, the expression for the velocity poten- 
tial takes the form of a power series in the transonic similarity 
parameter. On the basis of this form of the solution, an attempt 
is made to solve the nonlinear differential equation for transonic 
flow past the wavy wall. The analysis utilized exhibits clearly 
the difficulties inherent in nonlinear-flow problems. The inves- 
tigation has nevertheless been rigorously carried to the point 
where the question of the existence or nonexistence of a mixed 
potential flow free of discontinuities can be setiled by the behavior 
of a single power series in the transonie similarity parameter. 
The calculation of the coefficients of this dominant power series 
has been reduced to a routine computing problem by means of 
recursion formulas resulting from the solution of the differential 
equation and the boundary condition at the wall. One of the 
interesting results of the analysis is the rigorous statement that 
the transonic similarity parameter must be less than four-thirds. 


INTRODUCTION 


The present report considers in detail a notion first ex- 
pressed in reference 1, namely, that the Yon Kármán transonic 
similarity rule is implicitly contained in the potential-flow 
solution for high-subsonic flow past a prescribed body. The 
calculations are very much simplified by the choice of rec- 
tangular Cartesian coordinates as independent variables. For 
this reason the solid boundary chosen is a two-dimensional 
wavy wall of small amplitude extending to infinity in 
both the downstream and upstream directions. The problem 
is first treated by means of the Prandtl-Busemann iteration 
method for high-subsonic undisturbed speeds. On the basis of 
this solution an attempt is made to solve the nonlinear differen- 
tial equation for transonic flow (corresponding to the linear 
Tricomi equation in the hodograph plane) past a wavy wall. 
The analysis utilized exhibits clearly the difficulties inherent in 
nonlinear-flow problems. The investigation has nevertheless 
been rigorously carried to the point where the question of the 
existence or nonexistence of a mixed potential flow free of 
discontinuities can be settled by the behavior of a single power 
series in the transonic similarity parameter. The calculation 








of the coefficients of this dominant power series has been 
reduced to a routine computing problem by means of recur- 
sion formulas resulting from the solution of the differential 
equation and the boundary condition at the wavy wall. 

The author wishes to acknowledge the invaluable aid and 
advice of Dr. A. Busemann of the Langley Laboratory during 
the writing of this report, especially with regard to the final 
section “General Analysis.” 


CALCULATION OF HIGH-SUBSONIC FLOW PAST A, WAVY 
WALL BY MEANS OF PRANDTL-BUSEMANN 
ITERATION METHOD 


The fundamental nonlinear differential equation for the 
potential flow of a compressible fluid can be written as 


2 2 
Mau?) Oxx +a brr—2M o bxrýrxr =Q 
(1) 
where 
c? —1 5 
mia Viwe 2) 
and 
¢ velocity potential of flow 
A, Y rectangular Cartesian coordinates in plane of flow 
u,v fluid velocity components along X- and Y-axis, 
respectively 
U undisturbed stream velocity 
€ local speed of sound 
Co speed of sound in undisturbed fluid 
AM, Mach number of undisturbed stream (U/e,) 
Y ratio of specific heats at constant pressure and con- 


stant volume 
The quantities ġ, X, Y, u, and v are nondimensional with a 
characteristic length / as unit of length and the undisturbed 
stream velocity U as unit of velocity. The subseripts X and 
Y in equation (1) denote partial differentiation with respect 
to the designated variables. 

In order to obtain the Prandt!l-Busemann iteration equa- 
tions based on small perturbations of the undisturbed 
stream, the assumption is made that the velocity potential 
ọ can be expanded in the form 


=AX+Oit¢b2tOs+ o (3) 


t Supersedes NACA TN 2383, “On a Solution of the Nonlinear Differential Equation for Transonic Flow Past a Wave-Shaped Wall” by Cari Kapisan, 1951. 
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For the purpose of defining and controlling the iteration 
procedure, the function ¢,4, and its derivatives are re- 
garded as small compared with the preceding approxima- 
tion ¢@, and its derivatives. From equation (3) and the 
fact that for irrotational or potential flow u=—¢y and v=doy, 


u=1l+éixtéatésxet ... 


and 


w=dyy thr trt .. . 


When these expressions for u and v are introduced into 
equation (1), together with the corresponding expression 
for jea? given by equation (2), and the powers and prod- 
ucts of ¢, and their derivatives are grouped according to the 
assumptions of the small-perturbation method, the following 
iteration equations for the first three approximations result: 


duet diy =0 E 
buat omy = 2M N + o)budurt ónt] (5) 
burt dav 2A? { (1-0) Crab} didn) + 
[26° +e)—Idiedudiat 
Hof 280 +05 found? tHe Noida! + 





Dixy Pay sudan} (6) 
where z and y are new independent variables defined by the 
transformation 

r= X 
(7) 
y=8Y 
and 
e?=1—Af,,? 
yti Ma : 
sae at 


Equation (4) is a Laplace equation and equations (5) and 
(6) are Poisson equations where the right-hand sides contain 
only previously determined quantities. These equations 
have been treated recently in reference 2 where the particu- 
lar integrals of equations (5) and (6) are given in real form. 
These particular integrals have been utilized in obtaining 
the flow over the wavy wall. Thus, the equation of the 
infinitely long wave-shaped wall (fig. 1) is assumed to be 


Y=a cos aX (8) 
where 
a amplitude of wave 
see 
` 


A wave length 








y 











Y 


Ficure 1.—Wave-shaped wall. 


The reference element of length is conveniently chosen to 

1 A : . . 
be a (=>) The equation of the wavy wall in nondimen- 
sional form and in termsof the variables zand y then becomes 

y= ße cos & (9) 
where, if the thickness coefficient of the wavy wall is defined 
a 
as t= ———, then e=t. 
Qe 

The expression for the velocity potential ¢, obtained with 
the aid of the particular integrals given in reference 2, 
including terms of the order é and satisfying the boundary 
conditions at infinity and at the solid wall to the same order 
is as follows: 


€ 


paat$er sin eg (F) (2eMay+ tet 
(2—0)B4} e-™ sin 224 (5) « 84M. [24+20— 


Llo2—o (22+13¢) 84} e+ M,?[o (6+50)+ 
(8+100+30%) BJ e-* +40 (2+0) Matye) sin a+ 


“5 9) { (27+ 4Al,?(—25+11¢0+40?—96%)] e- 


6 AL,,2 (3 +20+38% ye-#+ 907M *y’e-*} sin 3st... 
(10) 


It is of interest to examine equation (10) when the as- 
sumptions of transonic-flow theory are introduced. These 
assumptions are essentially that the undisturbed flow veloc- 
ity differs only slightly from the speed of sound, that the 
velocity component normal to the oncoming flow is small 
compared with the speed of sound, and that the velocity 
component in the direetion of the oncoming flow is of the 
order of the critical velocity c*. If the undisturbed stream 
is in the direction of the positive z-axis, then the velocity 
potential ¢, referred to the critical velocity e*, can be written 
as (see_reference 1) 


part -ig CMa) Sew) (11) 
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The second term on the right-hand side of this equation 
is the disturbance-velocity potential and implies that terms 
involving powers of 1—f,? higher than the first are to be 
neglected. The differential equation satisfied by the func- 
tion f(x,y) is obtained from equation (1) and takes the fol- 
lowing simplified nonlinear form: 


df af oF 


oy? Or Ox? (12) 








The boundary conditions to be fulfilled by f(x,y) are as 
follows: 


OF = 
TI. 
(at y= œ) 
oF _ 9 | (13) 
yY 
df as N si pams —o co 
7I sinz (at y=0, <r< | 


where /, the transonic similarity parameter, is 


2. rtl 

‘= My" 4) 
Thus, it is seen that f is a function of x, y, and the transonic 
similarity parameter $ only. In the limiting case, 1/,—1 
and e30, k retains its meaning as a transonic similarity 
parameter. 

Note that in the Prandtl-Busemann iteration procedure 
the order to which the boundary condition is applied at the 
surface of the profile is the same as the order to which the 
iteration has been taken. In the transonic case, however, 
the surface boundary condition is always applied at the 
axis y=0. In particular, this statement of the boundary 
condition for transoniec flow is a rigorous one in the limit as 
both the thickness coefficient and 1—J/,? simultaneously 
approach zero but with the same distribution of slope that 
belongs to the family of profiles being treated. A rather 
complete discussion of this point is contained in reference 3. 

Before application of the foregoing considerations to equa- 
tion (10), it should be noted that ¢ in that equation is re- 
ferred to the undisturbed stream velocity U. In order to 
introduce c* as the reference velocity, both ¢ and the right- 


hand side of equation (10) must be multiplied by U/c*. 
#2 


: . 7 : cet 
Now, from equation (2), with e=c* and wae it 


follows that 





1 Ma? 
M= (15) 
147 Ma 


where 


s U 








If terms of only the first power of 1— M? are retained, it 
follows from equation (15) that 


Mtsi1-—5 G17 5%.., (16) 


Then, multiplying the right-hand side of equation (10) by 

this expression for M*, replacing e/8 by at (—AL 7k 
T 

from equation (14), and neglecting all terms containing 

1—Ma to a power higher than the first yields the following 

expression for ¢: 


i R 1l f 
$=1+_ I (ies AF) f- eteta r) e~¥ sin r+ 


1 ; a 
T6 (1+2y) k’e-* sin 224| —355 (6+4y) sin r+ 
1 





(4+12y+9y?) sin az | keh tee (17) 


bo 


88 


Comparison of this equation with equation (11) shows that 
the result obtained by means of the Prandtl-Busemann 
iteration method contains implicitly the form of the solution 
required by transonic-flow theory. Moreover, the expres- 
sion for the function f(x,y), namely, 


fe, =—2+(k-+ 555 r) e™ sin «+ 


1 


(1+2y)k?e-™ sin 2z+| — 555 (5+4y) sin z+ 


= 
16 
1 


ggg (4+ 12yt 9y? sin 3r] ke-wy n.. (18) 


satisfies the nonlinear differential equation (12) and the 
boundary conditions, equation (13), to the order $°. This 
fact suggests a solution for the flow over the wavy wall, in 
the neighborhood. of Mach number unity, obtained directly 
from equation (12). 


SOLUTION OF EQUATION FOR TRANSONIC FLOW PAST 
A WAVY WALL 


Equation (18) suggests the following form for a solution 
of the nonlinear differential equation (12) for transonic flow 
subject to the boundary conditions stated in equation (13): 


f(z, y= —t+>) Jasin nz (19) 
where the/,’s are functions of y only. This form for f(x, y) 


is substituted into equation (12) and repeated use made of 
the following identity: 
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When the coefficients of the separate harmonic terms sin nz 
are placed equal to zero, the following system of nonlinear 
ordinary differential equations for f, results: 


fa! nfr = -} n > mn—mM)fmfrom— 


l e2] 
2 n Day (m+ 1)(m +n Din iil nent 
m= 


(n=1,2,... @)} (20) 

Note that equation (20), with the right-hand side placed 
equal to zero, has a solution of the form e7™. An iteration 
procedure is then set up in such a manner that the highest 
power of k and of e™” appearing in the functions f, is equal to 
the number designating the iteration step. Thus, the fol- 
lowing sets of iteration equations result: 


Step 1, 
fi’ —fi=0 (21) 
Step 2, 
fll Afra —5 f? (22) 
Step 3, 
A — h= fife } 
(23) 
J —9fs=— 3f f: i 
Step 4, 
fi! —4f=—} FP —3fifs 
(24) 
4 —16 fs —6 fi fof? 
Step 4, 
F — h= fifa 3 hfs 
R —Ifs=— 3 fifa bff (25) 
: Ss’ —25fs=—10fifi— 1 5 fof 
Step 6, 


fy —4f_= -; Fe—3fi fa— Shots 
fe’ —16fs=— 6 fi fa—4f?—1 0S, fs (26) 


a 27 
fe! —36 fe=—1 5f fE f 


and so forth, 

The right-hand sides of these equations are known functions 
constructed from previously determined quantities in accord 
with the iteration procedure adopted. 

These equations are of the second order, nonhomogeneous 
type with constant coefficients and are readily integrated. 
The resulting expressions for fi, fo, Js, Ja fs, and fs, with the 
boundary conditions (13) taken into account, are as follows: 


24 fof. 


1861 Nc Ger 2765 a, 
i= =(14+55; K Eeoae i oto naset $ 
a eee 65 
G 64 txa) y |e i -a0 tja y+ 
9 Hie i °) 355 —5y¥ 1 rå 
32x256 ” Fason! Be ey (27) 
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yt) He, (28) 





47 34 
192X256 y Isasa! 
1 — 23603 


EEN Eza rs ae 259 k? 
fe=| ab 97056" +(satyascone t v+ 


Ly 2 a ET 155 
(z+ grag) ‘|i "3X256 32256 YT 











T17 ag b +) 15-5 
3952256 Y +23 a (29) 
390547 _,, 32947 e) 
t= otraa et 384" a y+ 


3 617 2 i) | wt AY N 
125 1288x256" ‘yy +t IZo ene 





35251 2465 455 

JOX256 1576X256 ” 96x256” T E 
t19 s 25 S ‘) Ten an 

73256 ” tas yf) ert (30) 
=t +765 vt Y air yt 

5,5: 
sascae v*) He EEs (31) 
91 91 23, s 
J=( 541x2045 "72x256 YT 3048 Y tgar er 
15 
2048 ” ttig vt) e+.. (32) 


Note that the functions fı, fe, and fs include the terms of 
equation, (18), obtained from equation (10) by allowing the 
Mach number to approach unity. 

Equations (27) to (32) may be considered to be essentially 
the nonlinear solution for the flow past a wavy wall of small 
amplitude for stream Mach numbers in the neighborhood of 
unity in the form of a power series in the transonic similarity 
parameter k. Moreover, this solution is identical with the 
one obtained by means of the Prandtl-Busemann iteration 
equations when Mach number unity is approached. 


CALCULATION OF LOCAL MACH NUMBER, CRITICAL 
SIMILARITY PARAMETER, AND 
PRESSURE COEFFICIENT 


GENERAL FORMULAS 


From equation (11), when all terms containing 1—A/.? 
to a power higher than the first are neglected, the expression 
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for the fluid velocity referred to the critical speed of sound 


$ 


c* is 


e=1t—2 0M. 52 (33) 


Now, from equation (2), the relation between q and the 
local Mach number Af is 


Lag r 


5 
= (34) 
1+% lyp 





2 


c= 





In the transonic approximation, the difference of any Mach 
number in the field of fow from unity is considered a small 
quantity. If terms of only the first power of 1—Af are 
retained, equation (34) yields 





G=1+ 7, AF-D)... (35) 
Hence, from equation (33), 
1— M= Me (36) 


This equation provides a means of calculating the critical 
value of the transonic similarity parameter; that is, the 
value of the parameter k for which Af=1 at the point of 
maximum velocity on the boundary. 


CALCULATION OF THE CRITICAL VALUE OF k 


For the family of wavy walls of small amplitude (including 
the limiting case of vanishingly small amplitude) with Jf=1 
at z=0 and y=0, equation (36) yields the following relation 
for the determination of the critical value of k: 


af 2 
=e) =0 (37) 


¥v=0 

By means of equation (19) together with the expression 

for fı to fs given by equations (27) to (32), equation (37) 

yields the following power series, exact to seven terms, for 

the determination of the critical value of & for the family 
of wavy profiles: 








sap es a3 337 a, 4043 psa 
SETAE 256X576" ! 

359381 3, L 

a70x250" Te 5l (38) 


The procedure adopted in order to estimate the critical 
value of & is as follows: From equation (38) the value of 
k can be found for 2, 3, 4, 5, 6, and 7 terms. These values 
of & are, respectively, 1, 0.8990, 0.8644, 0.8504, 0.8424, and 
0.8377. The last two values indicate the approach to the 


asymptotic value of k, that is, the value of k when the number f 


of terms in equation (38) is infinite. If the values of k 
approach smoothly to the asymptote, the estimated critical 
value £=0.8377 is very nearly correct. 


272483 —54——16 





Suppose now that both sides of equation (36) are divided 
by [(v+1)e?*. Then, since 


equation (36) can be written as 


1-3? 
C+D 


The right-hand side of equation (39) is a function of x, y, 
and the parameter k only and is characteristic of the entire 
family of boundary profiles. For the family of wavy walls 
with the critical value assigned to &, equation (39) evaluated 
at the wall becomes 


1—3ar 
(Cy + 1) el" 


of 
2-2/3 L2 
a (39) 


=1.1253(1—0.8536 cos x —0.0989 cos 24 — 
0.0299 cos 3xr—0.0122 cos 4r—0.0038 cos 5x— 
0.0017 cos 6+) (40) 


Table I lists the values of ii for values of x be- 


tween —r and r and figure 2 shows the corresponding 
curve. 


CALCULATION OF THE PRESSURE COEFFICIENT 


Bernoulli’s theorem for a compressible fluid assumes the 
following differential form along a streamline: 


2 

Up—p.)-+5 Ud (f1 )=0 (41) 
where 
p pressure in fluid 
p density of fluid 
q speed of fluid 
and the subseript œ denotes the quantity in the undisturbed 
fluid. Now, from equation (2) the relation between the 
nondimensional speed g/U and the local Mach number M is 
given by 


@ eat = Ar 
me 1475+ Lap Ae? 





(42) 


In accordance with the assumption of transonic-flow theory 
that all Mach numbers in the flow differ only slightly from 
unity, equation (42) becomes 


2 
ii ae [1 


=f 2— 
U? yl Al 


(i—AL +... (43) 
where powers of 1—3 and 1—Af higher than the first 
have been neglected. Then with the pressure coefficient 
defined as 
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Ficure 2,—Distribution of 3 and 
equation (41} takes the following form for transonic flow: 


(1—Af,,%)] (44) 





(ape a ear 
=e 


After integration, 


2 


a Haia 
o=F7_ l AP—i-—M,,”)] 


CiM 
or, with the aid of equation (36), 

l 2 PE Of 

? AAEE AES — he 2 > 


The right-hand side of this equation can be considered to be 
the first term in a power-series development of C,, Me in 





1—Af,”. In particular, when the local Mach number first 
attains unity, then oF at z=0, y=0 and equation (45) 
becomes F : 
2 ; 
Catam EI (1 ~Me’) (46) 


a result valid in the transonic range only. Again, if the 
thickness coefficient approaches zero as Af,—-1 and M=1 
at z=0, y=0, then equation (46) shows that the slope 



































Cy, ,, at the surface for the family of wave-shaped walls. 


Ja- M ®t the critical value Mu =Ma >l is a constant 


2 
yl 
treated. This result is valid whether the approach to Mach 
number unity is made from the subsonic or supersonic region. 





» independent of the particular family of profiles 


1/3 

If both sides of equation (45) are multiplied by ot, 
then 

1 173 ò z 

ae? Cie ee (1432) (47) 


where the right-hand side depends only on z, y, and the 
parameter & and is characteristic of the family of boundary 
profiles treated. For the family of wavy walls with the 
critical value of k chosen, equation (47) takes the following 
form at the solid surface (y=0): 


1/3 
ee Cp, M., =—2.2507(0.8536 cos 2+-0.0989 cos 22-4 


0.0299 cos 32-+0.0122 cos 42-+ 
0.0038 cos 52+0.0017 cos 62) (48) 


1/3 
Table I lists the values of arp" Cp, Mer for values of x 


between —7 and r and figure 2 shows the corresponding curve. 
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— Afz a v3 
TABLE I—VALUES OF jfr- wa oe AND — Crue, AT 


THE SURFACE FOR THE FAMILY OF WAVE-SHAPED WALLS 

















i | 
H 3 : H 
i og 1—4 EDS 
! (deg) [y+ 1) ere a 
—150 1.998929 1.743190 
175 1. 994505 1. 738342 
—170 1. 987188 1. 723708 
—160 1. 957430 1. 664191 
150 1. 906695 1, 562722 
—140 1, 834890 1.419111 
~130 L 743433 1. 236197 
—120 1. 634668 1. 018664 
i —110 1. 510815 770962 
100 i 1. 373602 498536 
-95 | 1, 300524 . 350379 
—90 1. 224843 199018 
—85 1. 146916 043163 
50 1. 067146 —. 116375 
-70 - 903656 —. 443357 
—60 -737110 —. TTH 
—50 - 569667 L1. 111335 
—40 -404761 —1. 441145 
| —30 < 250: —1. 750223 
—20 . 119935 —2. 010799; 
i —10 031487 —2. 187694 
| -5 i -007977 —2. 234715 
i a i 000007 —2. 250654 
i 5 -007977 —2 234715 | 
! 10 031487 —2. 187694 i 
20 . 119935 —2 010799 
30 1250223 —1. 750223 i 
i 40 | -404761 —1. 441145 
50 . 569667 —1. 111335 ; 
0 . 737110 —. 776449 : 
70 . 903656 —. 443357 
£0 1.067146 —. 116375 
85 1. 146916 . 043163 | 
90 1, 224843 . 199018 
95 1. 300524 -350379 
} 100 1.373602 . 498336 | 
i 110 1. 510815 < 770962 
| 120 1. 634666 1. 018664 
130 1. 743433 1. 236197 
i 140 1, 834890 3. 419111 
150 i 1. 906695 1. 526722 
180 1. 957430 1. 664191 
I7 i 1. 987188 1. 723708 a 
130 l 1. 996929 1, 743190 








GENERAL ANALYSIS 


An examination of the expressions for f, to fe given by 
equations (27) to (32) shows that the general form of fa is 


B 


fa=2 e 


p=0 


g= 
~ (2p +n)y yey Ate eee 
q=0 p ' 
(n=1,2,... œ) (49) 


where, if p=0, the upper limit of q is n—1, and, if p0, the 
upper limit of q is 2p-++-n—2. The four-labeled coefficients 
cit? are real numbers calculated from recursion formulas 
obtained from the system of differential equations (20) and 
the boundary condition at the surface of the wavy wall. 
The boundary condition at y= œ is automatically satisfied 
by the form of f,; whereas the boundary condition at the wall 


takes the form 
: =Í) 5 
(fa lr=0= 7 = (50) 


Inserting the expression for fa given by equation (49) into 
equation (50) yields immediately the following results: 


AM 
and 
r r 
a (2p+n) d3? = mm Ae (n=1,2, E GES a œ) (51) 
p=l = 


where, if n=1, the lower limit of p on the right-hand side is 
unity and, if n1, the lower limit of p is zero. Also, if 





n=1, the upper limit r of p goes from 1 to œ and, ifuxl,r 
goes from 0 to œ. 

In order to find the recursion formulas for the coefficients 
At?, the expression for f, given by equation (49) is sub- , 
stituted into the system of differential equations (20). The 
calculation is facilitated by the introduction of the following 


notations: 


anaa e 
A= Dy yo Ape (52) 
g=0 r=p 
and 
A=" => A® wes (5 3) 
g=0 


where the quantities A*,* arise from the multiplication of 
the two infinite series >} Ate?” and >) Aze, Note 
p=0 p=0 


that the quantities 4%” are symmetric with respect to the 
upper labels m and n. Then 


Famke Ax eu? (54) 


When this expression for f, is substituted into equation (20) 
and repeated use is made of the identity 


the exponential terms in y can be eliminated, and the follow- 
ing recursion formulas result: 


n=l 
—2n(A3)' += nS) mn= mAy 


(p=0; n=2,3,... 0) (55) 
and 
4plp tnd —2(2p +n + y= 
jn = m(n—m)A™*, eee = mint m)k?™ Amt tm 
(p=1,2, . eyn=—1,2,... 0) (56) 


Note that these recursion formulas still contain powers of y 
and k. Both y and $ can be eliminated for given values of p 
and recursion formulas containing only the coefficients 
Az? thus obtained. For example, consider equation (55) for 
which p=0. By repeated use of the identity 


n a n n 
Dy Dy Bas= 2, 2, Ba,p 
a=0 £=0 B=0 a=8 
the following recursion formula is obtained: 


—2nl¢ +A sie t@+ D@t2e 7A j= 


q n—2-i 
zi >i mt Dm- m— 1) = Am eae ped 
4 t=0 m=g— 
(n=2, 3, o;g=0,1,2,...m—2); 
r=0,1,2, œ) (57) 
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where 


A number of interesting relations can be easily obtained 
from this recursion formula. Thus, for g=n—2 and r=0, 


20—14, $F D asm) A R AMES 
m=0 
Then 

420l so _ 25 

di g A o 48X128 
i 9 
1 

A3 o= 


From these numerical values the general form is found to be 


-9 


n” 2 


A,” =p 
Similarly, for g=n—2 and r=}, 
A,” Lately 


In a corresponding manner, more complicated expressions 
can be obtained for A,%, ° 0 whan r=2, 3, : 

Note that in the expression for the distumbancesnouential 
function f(x,y), infinite series of the type 


Syyte Mh"? sin nz A,r, 8 
=1 
The ratio of the (rn-+-1)st and nth terms is 
sin (m+1)z Arti 
sinne  A,%,° 
In order for the series to converge, the limit of this ratio as 
n—>æ must be less than unity. Thus, if the maximum 


value of ye"(=2) is inserted and #0, Cauchy’s ratio 


occur. 


kye™ 





test yields 
ke É 
Saol or Å< 4 (59) 
Other infinite series oceur in the expression for f(e,y) which 
diverge for values of the transonic similarity parmeter + 
considerably less than 4. Thus, consider the recursion 
equation (57) and take g=0 and r=0. Then 
2nA? 04255? 50 
na 
-3n 5 (M+ 1)(n—m DAMS 3 ANH 
m= 





and 


ot 2( AG 0)? 
=2 A3 0+ 10450+1543 3 
124} 3=2 4$ 9-1543 +2443 ASIHE (4D 
Also from equation (57), 
443 0= Aid 
84i 2=34 
a 
843 0 =243 0 +8438 AG 


390+15439 AZ 0+104i3 
ot GATS Aoo t 8A n A oto] o 


and 


3043 0=12 A1 H15 A30 A30 t10 Aio t15 Ao ATS 


124$ 9=4 Al 3+9 433 A334344 A 8+5 ASS + 


s ALBEE (AID? 


20 
244 20 AT 


o deot 15 A30 +27 Ad dio 


e.. sa 


The following relations are obtained from the supplementary 
equations (51): 


From these relations, together with the ones listed in equa- 
tion (58), 





1 7 
At =g Aš o=768 
i 7X13 
A= 3y At =S E 
7 
A} °—365 eee ee 
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and lead to the general rule 


no. on—5} 
at Sema ad (n=2, 3, wee œ) (60) 
From equation (51), nA §=<7 8; therefore, 
„o [3n—5} 
A =a T (n=2, 3, oe * oo) (61) 


where, by definition, 


(3n—5}=1X4X7X10X13X ... X(8n—5) 


In the expression for the local Mach number in the field of 
flow given by equation (36), the expression 0f/Oxr occurs. 
This expression, with z=0, contains infinite series of the 
type 


F=3) nk*e~™ Ase (62) 
n=l 
and 
G= nk"e™ Ars (63) 
n=2 


If the maximum value of e~*#(=1) and the expressions given 
by equations (60) and (61) are inserted, the Cauchy ratio 


test shows that both series converge for kK In particular, 


the series expression for F evaluated at the surface (y=0) 
can be expressed in closed form. Thus, 


Pepi ett pe 2 (1-4 (64) 


n= ni4r- i 


The graph of F against k is a semicubical parabola with the 


cusp point at p= and F=2. With the restriction that the 


transonic similarity parameter k be positive and that one 
and only one value of k correspond to a given value of F, 
the permissible values of k and F are confined to the part 
of the parabola lying between the origin (0,0) and the apex 
(4/3,2). The power-series expression for Ẹ@ evaluated at 
the surface can also be expressed in closed form; namely, 








This expression, together with the one for F, shows clearly 
that the parameter k cannot be equal to but must be less than 
four-thirds. 

A close examination of the recursion formulas (55) and 
(56) discloses the important fact that each one of the mani- 
fold of power series in $ that appear in the functions fa can 








be expressed in terms of the members of a single dominant 
set of power series. This dominant set consists of one power 
series in k from each f,, namely, the one multiplied by e7" 
only. According to equation (49), the members of the 

dominant set of power series are given by 
S, =9 Ago kets (n=1,2,... œ) (66) 

T=Q 

Several examples are now given to illustrate this impor- 
tant observation. Consider the series that belong to the set 


a 
2, dpp bet 


r=0 


(n=2, 3,... œ) 


From the recursion formulas, 





r=0 
= i 1 1 
Dy Ate ke ier 2S a+ Ž SS+ Sr 
Consider now the set 
Dy di? krata (n=3, 4, œ) 
Then 
= i 1 
x Ag? het S3 
z At? kYti=— 1 setiges 
2 ir =Jog V TJ >92 
Consider the set 
St Agl kate (n=1,2,... ©) 
r=] $ 
From the recursion formula (56), 
2 Al} krt l=—525S, 3 -5 S:S: 
Sy Ag? ta Si StS SS 
Z Or 12x256 1 imag Pies 
= ; 163 
ASETTI m 2 
2 dirk 2569" -ip 5 È 355 5! $S — 
5 S S238, S: 


s.r’ 
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Consider the set— 


(n=1,2,... &) 


T=] 
Then 
5> Ai} prL g 3 
— lr 64 1 


© 3 } 
pa At) prt? OG tT Sg, 
2y Aik 356 © ye 





D AN ere TE SY PE SiS SSF SS? 

Other examples are 

= ie Ss 

2 AN Haig St— yga 5*9 an 2 Gz SiS 

Poy Ade RP E a Si P SiS gg 5S4 l 
5 818.4 SaS; 


The recognition of the existence of a dominant set of power 
series in & represents a major reduction in the complexity of 
the present problem. Thus, consider the array of infinite 
power series contained in equation (66): 


Si ARSE+ALS REALS ESL ALS RTH 
Sp AR KHAR? RT ARS RL ARS ASL 
S= Ab 0 kART P+ AGS k HAS ket 
SASS AAS! RL ASS ES ASS EH 
Spa A59 MASS MLAS OL ARLES 


(67) 


An examination of equations (27) to (32) shows that the 
coefficients Aj? of the series S, appear to be positive and 
monotonically decreasing. The series formed from the first 
column on the right-hand side of equation (67) therefore 
dominates the series formed from succeeding columns. 


Moreover, the first-column series 2433 " has a radius of 


convergence kas (see equation (61}). The radii of con- 


vergence of the other columnar series therefore are either 
equal to or greater than four-thirds. Similarly, an examina- 
tion of the series S, shows that the coefficients in each column 
on the right-hand side of equation (67) also seem to form 
positive and monotonically decreasing sequenees. This be- 
havior means that S, is the dominant series of the set S, and, 








in fact, of the aggregate of power series in / in the expression 
for the disturbance potential f(z, y). 
Consider now the series consisting of the first terms of the 


odd-labeled series S,, S3, Ss, . . . , that is, 
TAa ķ2n+ti (68) 


According to the theory of power series (and it can be easily 
verified), the radius of convergence of this power series is 


still =$. 


Now, a comparison of corresponding terms with 


the dominant series S, shows that 


and 
1861 


E ean T 
433256 At TE a5E 


Thus, if in general Aj ?>A?,+13 then the radius of conver- 
gence of the dominant series S, can be less than the radius 
of convergence of the comparison series given by expression 
(68) and therefore may conccivably be equal to the critical 
value k,,=0.8377. Moreover, it would then follow that the 
original Prandtl-Busemann small-perturbation method is 
valid for purely subsonic flows only. This conclusion would 
not invalidate other approaches to the transonic-flow problem 
(reference 4). 

Unfortunately, the coefficients 4} 2 do not conform to any 
apparent or superficial law, but perhaps a careful study of 
the recursion formulas (55) and (56) and the supplementary 
or boundary relations (51) will yield a rigorous proof of the 
foregoing statements. Otherwise, it remains to calculate a 
reasonable number of the coeficients A}?. For this pur- 
pose the development of complete recursion formulas similar 
to equation (57) for the required values of p is worth while. 
Thus, for p=1, the recursion formula is 


Ant A, 2-20-42) (q+) An hs 


+(¢+1)(¢+2) Jë Pi T Ata = 
1 ae 


=0 m 


ee m—1) 


| 


r 
PAA ek mam) Ag hae 8 |- 
$= 


5 mnt) dy Aste pete 
(n=1,2,...0; ¢=0,1,2,...; r=0,1,2,... 0) (69) 
where 
tin 0 (g=n—1 or n) 
Se (¢#n—1 orn) 


Finally, it may be of interest to give the general formula 
for the Mach number distribution at the surface of the wavy 
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wall. Thus, by means of equations (36) and (49), 


AL—1 


mela m 
2ml 55k" >) (n—2m) cos (n—2m)e DAV 
1—Af,? m=0 p=0 


n=l 
n : 
where [z] denotes the integral part of n/2. 
LANGLEY AERONAUTICAL LABORATORY, 


NATIONAL Apvisory COMMITTEE FOR AERONAUTICS, 
LANGLEY FIELD, Va., April 26, 1951. 


(70) 
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